arXiv: 1504.01141v3 [math.CA] 29 Oct 2015 


Generalized Chebyshev polynomials of the second 

kind 

Mohammad A. ALQUDAH 
October 30, 2015 


Abstract 

We characterize the generalized Chebyshev polynomials of the second 
kind (Chebyshev-II), then we provide a closed form of the generalized 
Chebyshev-II polynomials using Bernstein basis. These polynomials can 
be nsed to describe the approximation of continuons functions by Cheby¬ 
shev interpolation and Chebyshev series and how to compute efficiently 
such approximations. We conclude the paper with some results concerning 
integrals of the generalized Chebyshev-ff and Bernstein polynomials. 
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1 Introduction, background and motivation 

Orthogonal polynomials are very important and serve to approximate other 
functions, where the most commonly used orthogonal polynomials are the clas¬ 
sical orthogonal polynomials. The field of classical orthogonal polynomials de¬ 
veloped in the late 19th century from a study of continued fractions by P.L. 
Chebyshev. 

We have seen the significance of orthogonal polynomials, particularly in the 
solution of systems of linear equations and in the least-squares approximations. 
Meanwhile, polynomials can be represented in many different bases such as the 
monomial powers, Chebyshev, Bernstein, and Hermite bases forms. Every form 
of polynomial basis has its advantages, and sometimes disadvantages. Many 
problems can be solved and many difficulties can be removed by suitable choice 
of basis. 

In this paper we characterize the generalized Chebyshev polynomials of the 
second kind (Chebyshev-II) These polynomials can be used de¬ 

scribe the approximation of continuous functions by Chebyshev interpolation 
and Chebyshev series and how to compute efficiently such approximations. 
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1.1 Bernstein polynomials 

We recall a very concise overview of well-known results on Bernstein polynomi¬ 
als, followed by a brief summary of important properties. 

Definition 1.1. The n -|- 1 Bernstein polynomials of degree n, x G 

[0,1], fc = 0,1,..., n, are defined by: 


Bi(x) = 




— k 


fc = 0,1,..., n 
else 


( 1 . 1 ) 


where the binomial coefficients 


n\ 


k\{n — k)V 


k = 0,... ,n. 


The Bernstein polynomials have been studied thoroughly and there are a fair 
amount of literature on these polynomials, they are known for their geometric 
properties mu], and the Bernstein basis form is known to be optimally stable. 
They are all non-negative, B]^{x) > 0, x G [0,1], satisfy symmetry relation 
B^{x) = — x), and the product of two Bernstein polynomials is also a 

Bernstein polynomial which is given by 


B^{x)B^{x) 



-r>n-\-m 

^i+j 


(x). 


The Bernstein polynomials of degree n can be defined by combining two 
Bernstein polynomials of degree n — 1, where the kth nth-degree Bernstein 
polynomial can be written by the known recurrence relation as 


-Sfc (a;) = (1 - ^{x) + xBf:_l{x), fc = 0,..., n; n > 1 

where Bq(x) = 0 and B^(x) = 0 for k < 0 or k > n. Moreover, it is possible to 
write each Bernstein polynomials of degree r where r < n in terms of Bernstein 
polynomials of degree n using the following degree elevation [3]: 


Bl{x) = 


n—r-\-k 

E 

i—k 


OiTl) 

C) 


k = 0,l,...,T. 


( 1 . 2 ) 


In addition, the Bernstein polynomials can be differentiated and integrated eas¬ 
ily as 

■^B]^{x)=n[B^Zl{x)-B^-^{x)], n>l, 


and 


BJt{x)dx = 


n -I- 1 


, fc = 0,1,..., n. 


(1.3) 


These analytic and geometric properties of Bernstein polynomials with the 
advent of computer graphics made Bernstein polynomials important in the 
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form of Bezier curves and Bezier surfaces in Computer Aided Geometric De¬ 
sign (CAGD). The Bernstein polynomials are the standard basis for the Bezier 
representations of curves and surfaces in GAGD. 

However, the Bernstein polynomials are not orthogonal and could not be 
used effectively in the least-squares approximation m- Since then a theory 
of approximation has been developed and many approximation methods have 
been introduced and analyzed. The method of least squares approximation 
accompanied by orthogonal polynomials is one of these approximation methods. 


1.2 Least-square approximation 

The idea of least squares can be applied to approximating a given function by 
a weighted sum of other functions. The best approximation can be defined 
as that which minimizes the difference between the original function and the 
approximation; for a least-squares approach the quality of the approximation is 
measured in terms of the squared differences between the two. The following will 
briefly refresh our background information to enable us to combine the superior 
least square performance of the generalized Ghebyshev-II polynomials with the 
geometric insight of the Bernstein form. 

Definition 1.2. For a function f{x), continuous on [0,1] the least square ap¬ 
proximation requires finding a polynomial (Least-Squares Polynomial) 

Pn{x) = aoipo{x) -f ai(pi{x) H-h a„(pn{x) 


that minimize the error 


E{ao,ai,...,an) = / [f {x) - Pnix)]"^ dx. 


For minimization, the partial derivatives must satisfy 


dE 

doi 


= 0, i = 0,..., n. 


These conditions will give rise to a system of n -I- 1 normal equations in n -|- 1 
unknowns: oq, oi,..., a„. Solution of these equations will yield the unknowns: 
ao,ai,... ,a„ of the least-squares polynomial Pn(x). It is important to choose 
a suitable basis, for example by choosing y}i(x) = x^, the matrix coefficients of 
the normal equations given as 


^0 

which is Hilbert matrix that has round-off error difficulties and notoriously ill- 
conditioned for even modest values of n. 

However, the computations can be made efficient by using orthogonal poly¬ 
nomials. Ghoosing {ipQ{x),ipi{x),... ,y}n{x)} to be orthogonal simplifies the 
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least-squares approximation procedures. The matrix of the normal equations 
will be diagonal which simplifies calculations and gives a compact closed form 
for Oi, f = 0,1,..., n. 

Moreover, knowing Pn(x) will be sufficient to compute a„+i to get Pn+i{x). 
See [10] for more details on the least squares approximations. 


1.3 Factorial minus half 

We present some results concerning factorials, double factorials and some com¬ 
binatorial identities. The double factorial of an integer n is given by 

(2n — 1)!! = (2n — l)(2n — 3)(2n — 5)... (3)(1) if n is odd 
n!! = (n)(n — 2)(n — 4)... (4)(2) if n is even, 

where 0!! = (—1)!! = 1. 

From (II.4p . we can derive the following relation for factorials 

if n is even 

if n is odd • (^■^) 

From dm we obtain 

(2n)!! = [2(n)][2(n - 1)]... [2.1] = 2"n!, (1.6) 



and 


(2n)! = [(2n-l)(2n-3)... (1)] ([2(n)][2(n - l)][2(n - 2)]... [2(1)]) = (2n-l)!!2"n 

(1.7) 

By combining dmi) and dni), we get 


2n 

n 


■ 


In addition, using dm and with some simplifications we obtain 

ill) (2n-l)!! 

(^) (2fc-l)!!(2n-2/c-l)!!- 


( 1 . 8 ) 


(1.9) 


1.4 Univariate Chebyshev-II polynomials 

Let us hrst consider a definition and some properties of the univariate Chebyshev 
polynomials of the second kind. 

Definition 1.3. (Chebyshev polynomials of the second kind Un{x)). The Cheby¬ 
shev polynomial of the second kind of order n is defined as follows: 


= „ = (1,10) 


5277 COS 


Kx)] 
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From this definition the following property is evident: 


Un{x) 


sin(n +1)0 
sin9 


X = cos 9. 


( 1 . 11 ) 


The Chebyshev polynomials are special cases of Jacobi polynomials {x), 
and related as ^ 

17„(x) = (n+l)(^”+2^ (1.12) 

Authors are not uniform in orthogonal polynomials notations, and for the 
convenience we recall the following explicit expressions for univariate Chebyshev- 
II polynomials of degree n in x, see Szego m-- 




(2n+l)!! 


A;=0 


1 
2 

n — k 


X -\-l 


i—k 


X —1 


(1.13) 


which it can be transformed in terms of Bernstein basis on a; e [0,1], 




(2n+ 1)!! 




o 


(1,14) 


Also, it may be represented in terms of Gauss hypergeometric series as follows 

0 : 

Ur,{x) := (n + 1)2Fi | . (1.15) 

Although the Pochhammer symbol is more appropriate, but the combinatorial 
notation gives more compact and clear formulas, these have also been used by 
Szego m- 


1.4.1 Properties of the Chebyshev-II polynomials Un{x) 

The Chebyshev-II polynomials Un{x) of degree n are orthogonal polynomials, 
except for a constant factor, with respect to the weight function 

W{x) = -\/l — x^. 

In addition, Chebyshev-II polynomials satisfy the orthogonality relation [3] 

x^l-x)iUn{x)Um{x)dx = if m = u ■ 

The univariate classical orthogonal polynomials are traditional defined on [— 1,1], 
however, it is more convenient to use [0,1]. 
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2 Main results 


In this section, we characterize the generalized Chebyshev-II polynomials 
then we write them using Bernstein basis. Finally, we conclude the section with 
the explicit formula for the generalized Chebyshev-II polynomials using Bern¬ 
stein basis. 


2.1 Characterization 


Using the relation (11.121) and a construction similar to [51 [7] for M, N > 0, the 
generalized Chebyshev-II polynomials < ‘Wn^’^\x) \ can be written as 

t J n=0 

t Un(x)+MQ„(x)+NE„(x)+MNS„(x), n = 0,1,2,... 
z^(n + Ij! 

( 2 . 1 ) 


where for n = 1, 2, 3,... 


Qn{x^ — 
Rn{,x') — 


{2n+l)\\ 

3.2"n! 

(2n-h 1)!! 
3.2’^n! 


n{n + 2)Un{x) - -{x - l)Z)t7„(x) 


n{n + 2)Un{x) - -(x + l)DUn{x) 


and 


Su{x) = - l)i' 2){7„(x) - 3xDUn{x)]. 


( 2 . 2 ) 

(2.3) 

(2.4) 


If we use 

{x^ — l)D^Un{x) = n{n + 2)Un{x) — 3xDUn{x), n = 0,1, 2, 3,... (2.5) 

we easily find from dH that 


Sn{x) 


{n + 2)\{2n + ly.l ^ 
3^.2”n!(n — 1)! 


l)D^Un{x), 71=1,2,3,.... 


( 2 . 6 ) 


It is clear that Qo{x) = Ro{x) = So{x) = 0. 

We know that the generalized Chebyshev-II polynomials satisfy the symme¬ 
try relation [7], 


n = 0,1, 2,... 


(2.7) 


which implies that Qn{x) = (—l)"i?„(—x), 5'„(a:) = (—1)"«S'„(—x) for n = 

0 , 1 ,.... 

From (ESI) and (lOl) if follows that 


^ 71 ^ + 2)(2n + 1)!! rr 71 1 _ 1 o o 

Qn(l) — 3 2"(n-l)! n-1,2,3,... 


( 2 . 8 ) 
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(2.9) 




(n + 2)(2n + l)!! 
3.2"(n- 1)! 




n= 1,2,3,... 


Note that the representations (12.21) and (12.31) imply that for n = 1, 2, 3,..., 
we have 


I \ ^ - 4 .- 1 . n(n + l)(2n+l) in\ 

<ln = -s- (2-10) 


fc=0 


'2'=(fc + 1)! 


and 


D ! \ (2fc +1)!! rr /- ^ -j-n n(n+l)(2n + l) 

Rn{x) = ) , ^~,, Uk{x) with r„ =- - -. (2.11) 


fc=0 


'2fc(fc + 1)! 


We also can find from (12.41) that for n = 1,2,3,..., we have 


(2fc +1)!! rr -.U (n + 2)(n+l)2n2(n-l) 

S„{x) = } Sn = 


k=0 


'2fc(fc + 1)! 


( 2 . 12 ) 


Therefore, for M, N >0, the generalized Chebyshev-II polynomials 
are orthogonal on the interval [—1,1] with respect to the weight function 


OO 

n—0 


2 1 1 

— {1 — x)^ {1 + x)^ + Mi5(a: + 1) + A^(5(a: — 1). 

TT 

and can be written as 




(2n+ 1)!! 
2"(n + 1)! 


Un{x) + Afe 

fc=0 


(2fc+ 1)!! 

2'=(fc + l)! 


Uk{x), 


(2.13) 


(2.14) 


where 

Xk = Mqk + Nrk + MNsk- (2-15) 

The next theorem provides a closed form for generalized Chebyshev-II poly¬ 
nomial x) of degree r as a linear combination of the Bernstein polyno¬ 

mials i?[(a:), i = 0,1,..., r of degree r. Those results generalize [9] contributions 
to the univariate Chebyshev polynomials of the second kind. 

Theorem 2.1. ForM,N > 0, the generalized Chebyshev-II polynomials ‘^r^’^\x) 
of degree r have the following Bernstein representation: 


(2.16) 


2=0 




where Xk = Mqk Nrk MNsk and 

{2r + lf 


it- — 




2^''{2r -2i + l){2i + l) Q 


, i = 0,1,..., r. 
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where 


'^0,r — 


(2r + 1) f2r 
22 ’- 


The coefficients satisfy the recurrence relation 

_ (2r-2^ + 3) 

( 2 z + l) 


(2.17) 


Proof. To write a generalized Chebyshev-II polynomial [x) of degree r 

as a linear combination of the Bernstein polynomial basis Bl (a;), z = 0,1,..., r 
of degree r in explicit form, we begin with substituting (11.11) into (12.141) to get 




2=0 


fc=o ^ ^ i=o 


L-i) 


k— 


(2r + l)!!Y^ pr/ \ I \ (2^ + 1)!! i'|fc-z9 p^( \ 

' i=0 k=0 \ ^ 

(2.18) 


k=0 


j=0 


where 


79- — 

^i,r — 


(VV(Vf) . 


C) 


z = 0,1,... ,r. 


(2.19) 


Using (12.191) and applying (« + 5 )! = 2^^(27z + 1)!! with some simplifica¬ 
tions, we have 


r+h\fr + ^ 

r — i 


(2r-I- 1 ) (r - i)! (2r-|-1) (r - i)! 

(2r — 2z -I- l)(r — z)!(z — i)! (2z -|- l)z!(r — z — i)! 

2 *(2 r-hl)(r-l)!! 2’'-*(2r-h l)(2r - 1)!! 

“ 2’'(2r -2i + l)(r - z)!(z - 1)!! 2U!(2f -h l)(2r - 2i - 1)!! 

{2r + l) (2r-1)!! (2r-h 1) (2 r-l)!! 

“ 2^{2r -2i + l)(r - z)!z! (2i - 1)!! (29 -f 1) (2(r - i) - 1 )!!' 

Using the fact (2n)! = (2rz — l)!!2"7z! we get 


r + i 
^ 2 


(2r -I- 1)2 


22 ’-( 2 r -2i + l){2i + l)\r J \2iJ' 
For the recurrence relation, it is clear that for z = 1,..., r we have 




Thus, 




(--* + 1)0 ■ 
(2z + l)(2r +1)2(2;) (2Q 

22’'(2r -2i + l)(2f -h l)(2r -2i + 3)Q 


, z = l,...,r. 


( 2 . 20 ) 





















We conclude with an interesting integration property of the weighted gen¬ 
eralized Chebyshev-II polynomials with the Bernstein polynomials. To do this, 
we introduce the following definition. 


Definition 2.1. The Eulerian integral of the first kind is a function of two 
complex variables defined by 

B{x,y) = ( — uy~^du^ 5i(x), 5ft(?/) > 0. 

Jo 


Note that the Eulerian integral of the first kind is often called Beta integral. 
We observe that the beta integral is symmetric, a change of variables hy t = 1—u 
clearly illustrates this. 

Theorem 2.2. Let Bf{x) he the Bernstein polynomial of degree n and '^^ ^'^\x) 
be the generalized Chebyshev-II polynomial of degree i, then for i,r = 0,1,...,n 
we have 


x^{l-x)'^B'^{x)%^^'^\x)dx 


n\ {2i -I-1)!! 


rj 2*(z -h 1)! 


E(-i: 


i—k 




i + 


i — k 


B{r k A —— r — fc-l—) 




d=0 


j=0 


d i'i f d -\- 


J 


d-jj 


n B{r j -,n d - r - j -I -) 


( 2 . 21 ) 


where B(x,y) is the beta function. 
Proof. By using (I2.16L the integral 


/ = 



x^ (1 — x)^ Blf {x)‘^^^^’^\x)dx, 


can be simplified to 


1 = 


i2i + l)l\ C 

2*(f + 1)! Jo 


r-r+i 


(1 -x) 


n—r+7 








d =0 


(2d+l)!! 

2‘^{d+l)\Jo 




(1 -x) 


n-r+o 


1 / n 
r 


(D 

d (d+i)(d+i) 

j=o 


(2i -I- 1)!! /n 


2*(i -I-1)! \r 


E(-i) 


i — k 


k^O 


1+ /i + ^ 
i — k 


x^+f‘+i{l-x)^+^-^-'^+Ux 


Ea. 


(2d-hi)!! fn 


‘2'^(d-hl)! \rj ^ 
d=0 ^ \ 


E(-i) 


d-j 


d -\- f \ f d -\- 


d-j 


x^+J+2(l - a;)"+'^-’'-^+5dx 

( 2 . 22 ) 
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The integrals in the last equation are the Eulerian integral of the first kind, 
and can be written in term of Beta function as B{xi, yi) with xi = r + fc + |, 
j/i=n + i — r — A:+|, X 2 =r + j + |, and y 2 = n + d — r — j + ^. □ 
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